Dilatant strengthening as a mechanism for slow slip events

[1] The mechanics of slow slip events (SSE) in subduction zones remain unresolved. We suggest that SSE nucleate in areas of unstable friction under drained conditions, but as slip accelerates dilatancy reduces pore pressure p quenching instability. Competition between dilatant strengthening and thermal pressurization may control whether slip is slow or fast. We model SSE with 2‐D elasticity, rate‐state friction, and a dilatancy law where porosity evolves toward steady state ss over distance dc and ss = 0 + ln(v/v0); v is slip speed. We consider two diffusion models. Membrane diffusion (MD) is approximated by −(p − p)/tf where p and p are shear zone and remote pore pressure and tf is a characteristic diffusion time. Homogeneous diffusion (HD) accurately models fault‐normal flow with diffusivity chyd. For MD, linearized analysis defines a boundary E = 1 − a/b between slow and fast slip, where E ≡ f0 /bb(s − p), f0, a, and b are friction parameters and b is compressibility. When E < 1 − a/b slip accelerates to instability for sufficiently large faults, whereas for E > 1 − a/b slip speeds remain quasi‐static. For HD, Ep ≡ h/(b (s − p) ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi v1=chyddc p ) defines dilatancy efficiency, where h is shear zone thickness and v is plate velocity. SSE are favored by large h and low effective stress. The ratio Ep to thermal pressurization efficiency scales with 1/(s − p), so high p favors SSE, consistent with seismic observations. For Ep ∼ 10 transient slip rates, repeat times, average slip, and stress drops are comparable to field observations. Model updip propagation speeds are comparable to those observed along‐strike. Many simulations exhibit slow phases driven by steady downdip slip and faster phases that relax the accumulated stress. Model SSE accommodate only a fraction of plate motion; the remaining deficit must be accommodated during coseismic or postseismic slip.